arXiv:1503.03774vl [cond-mat.dis-nn] 12 Mar 2015 


Frontiers in Interdisciplinary Physics 


Research Article 
13 March 2015 


Record breaking bursts in a fiber bundie 
modei of creep rupture 

Zsuzsa Danku and Ferenc Kun"* 

Department of Theoretical Physics, University of Debrecen, P.O.Box: 5, Debrecen, 
H-4010, Hungary 

Correspondence*: 

Ferenc Kun 

Department of Theoretical Physics, University of Debrecen, P.O.Box: 5, Debrecen, 
H-4010, Hungary, ferenc.kun@science.unideb.hu 

Research Topic 


ABSTRACT 

We investigate the statistics of record breaking events in the time series of crackiing 
bursts in a fiber bundie modei of the creep rupture of heterogeneous materiais. In 
the modei fibers break due to two mechanisms: siowiy accumuiating damage triggers 
bursts of immediate breakings anaiogous to acoustic emissions in experiments. The 
rupture process acceierates such that the size of breaking avaianches increases whiie 
the waiting time between consecutive events decreases towards faiiure. Record events 
are defined as bursts which have a iarger size than aii previous events in the time series. 
We anaiyze the statistics of records focusing on the iimit of equai ioad sharing (ELS) of 
the modei and compare the resuits to the record statistics of sequences of independent 
identicaiiy distributed random variabies. Computer simuiations reveaied that the number 
of records grows with the iogarithm of the event number except for the ciose vicinity of 
macroscopic faiiure where an exponentiai dependence is evidenced. The two regimes 
can be attributed to the dominance of disorder with smaii burst sizes and to stress 
enhancements giving rise efficient triggering of extended bursts, respectiveiy. Both the 
size of records and the increments between consecutive record events are characterized 
by power iaw distributions with a common exponent 1.33 significantiy different from the 
usuai ELS burst size exponents of fiber bundies. The distribution of waiting times foiiows 
the same behavior, however, with two distinct exponents for iow and high ioads. Studying 
the evoiution of records we identify a ioad dependent characteristic scaie of the system 
which separates siow down and acceieration of record breaking as faiiure is approached. 

Keywords: Fracture, crackling noise, fiber bundle model, avalanche, record breaking statistics 


1 INTRODUCTION 


The fracture of heterogeneous materials procee ds in bursts genera ted by newly nucleating cracks or 
by intermittent propagation steps of crack fronts lATava et alJ (120061) . Measuring acoustic emissions of 
breaking bursts the fracture process can be decomposed into a time series o f crackling events , whic h 


contains valuab l e information about the microscopic dynamics of fracture iDeschanel et alj (120091) : 


Salminen et alJ (120021) : lRosti et alJ(l20Tol) . The analysis of crackling time series has usually been focused 


on the integrated statistics of events such as the probability distribution of the size (energy) and duration 
of bursts and of the waiting time between consecutive events. All these distributions are found to have 
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power law functional form wit h exponents having a hi g h degree of r obustn ess with respect to materials 


details ^schanel et alj(120091) : I Salminen et alJ l2002h : fRosti et alJ (I2010h . 


In applications materials are often subject to constant sub-critical loads which may lead to failure in 
a finite time. Such creep rupture processes are typical for components of engineering constructions and 
they also play a crucial role i n the emergenc e of natu ral catastrophe s such as land slides, stone a n d snow 
avalan c hes, an d earthquakes iNechad et alJ (I^05bll3): [Petriet al] (11994): iGarcimartin et alJ (Il997h 


Main 


(|20()0(): [Santucci et alJ (l2004linbes^anel et alj (l2006h : ISantucci et alj (l2006h : iKovfics et al. 


(120081) : iRamos et alj (12013l) . It is of high practical importance to understa nd how the creeping system 


approaches macrosco p ic failure using the data of acoustic moni toring IGarcimartin et alJ (119971) : 


Nechad et al.l (l2005bl) : iKovacs et al.l (120081) : iDeschanel et ^ (120091) . Computer simulations of discrete 
stochastic models of creep rupture are indispensable to analyze how the time series of crackling events 
evolves and to identify possible signatures of the imminent catastrophe. A unique feature of this evolution 
is that the rupture process is highly non-stationary, i.e. approaching failure larger and larger bursts are 
trigger ed while the proces s acce l erates indicated by the decreasing waiting time between consecutive 
events IGarcimartin et alj (Il997h : INechad et alJ (I2005bh. On the macro-sca le, the strain rate has been 


found to exhibit time-to-failure power law behavior [Nechad et alj (I2005bll3) : IDeschanel et alJ (I2006L 


20091) whic h is accompanied by the emergence of an Omori-type acceleration of the rate bursts on the 
micro-scale Danku and Kuii (l2013al) . 

In the present paper we investigate the evolution of the crackling time series of creep rupture by 
analyzing the statistics of record breaking (RB) bursts in a fiber bundle model (FBM) of creep failure. 
Records are bursts which have the larges t size since the b e ginning of the time series, hence, their 
behavior involves extreme value statistics iGalambo^ (119781) : lArnold et ^ (120111) . Motivated mainly 


by climate researc h i Redner and Petersenl (120061) and by the investi gation of earthquake time series 
Yoder et akl (|20 1 Ot) : Davids^ et aP (TiOOSh : [Shcherbakov et aO (120131) interesting analytical results have 
recently been obtained for the record breaking statistics of sequences of independent identically distributed 
(IID) random variables. The statistics of records ha s proven useful to i d entify t rends in tirn e serie s 
of measurements and to infer correlations of events I Shcherbakov et al^ (120131) : I Yoder et^ (120101) . 
Focusing on the limit of equal load sharing of our FBM we demonstrate that the record breaking statistics 
of crackling events provides novel insight into rupture phenomena. Comparing the outcomes of large scale 
computer simulations to the corresponding IID results on records we can identify regimes of the failure 
process dominated by the disorder of materials and by the enhanced triggering of breaking avalanches 
towards failure. The size distribution of records proved to have a power law form with a novel exponent of 
equal load sharing FBMs. Simulations revealed the emergence of a load dependent characteristic record 
rank which separates slow down and acceleration of record breaking when approaching failure. 


2 MATERIAL & METHODS 


To investigate the time series of breaking burst s we use a generic fiber bundle model (FBM) of 
creep ru pture introduced recently iKun et all (I2007LI2009L120081) : Halasz et alJ (120121) : Danku and KunI 
(l2013allbl) . We briefly summarize the main ingredients of the model construction emphasizing aspects 
most relevant for the present study. 


2.1 FIBER BUNDLE MODEL WITH TWO BREAKING MECHANISMS 

In the framework of the model the sample is discretized in terms of a bundle of parallel fibers which have 
a brittle response with identical Young modulus E. The bundle is subject to a constant external load ag 
below the fracture strength cxc of the system parallel to the fibers’ direction. Fibers are assumed to break 
due to two physical mechanisms: immediate breaking occurs when the local load cr* on fibers exceeds their 
fracture strength ah, i = 1,... ,N. Under a sub-critical load ao < ctc this breaking mechanism would 
lead to a partially failed configuration with an infinite lifetime. Time dependence is introduced in such a 
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way that those fibers, whieh remained intaet, undergo an aging proeess aeeumulating damage c(t). The 
damage meehanism represents the environmentally indueed slowly dev eloping aging of materi^s suc h 
as corrosion cracking and thermally or chemic ally activated degradation iKun et al] (12007112009L120081) : 


[Halasz et alJ (12012h : Danku and Kur 


similar to damage dynamics based models of rock fracture 


developed for instance in iTurcotte et alj (120031) . The rate of damage accumulation Ac* is assumed to have 
a power law dependence on the local load 


Aci = aajAt, 


( 1 ) 


where a is a constant and the exponent 7 controls the time scale of the aging process with 0 < 7 < + 00 . 
The total amount of damage Ci{t) accumulated up to time t is obtained by integrating over the entire 

loading history of fibers Ci{t) = a fg ai{t')'^dt'. Fibers can tolerate only a finite amount of damage so that 

when Ci{t) exceeds the local damage threshold the fiber breaks. Each breaking event is followed by 
a redistribution of load over the remaining intact fibers. Two limiting cases of load sharing are usually 
considered in FBMs: under equal load sharing (ELS) conditions all the intact fibers keep the same amount 
of load so that the load on a single fiber is ai = Nao/{N — i) after the breaking of i fibers. ELS 
realizes the mean field limit of FBMs where no stress heterogeneity can arise. In the opposite limit of 
localized load sharing (LLS) the load of broken fibers is equally redistributed solely over their intact 
nearest neighbors in the bundle. LLS leads to a high stress concentration around failed regions and it gives 
rise to the emergence of spatial correlations between consecutive breaking events. In the present study we 
consider only the ELS case where the homogeneous stress field hinders spatial correlation. In the model 


.th „th 


the quenched heterogeneity of materials is represented by the randomness of breaking thresholds cj 
i = 1,..., A. For simplicity, for both threshold values we assume uniform distributions between zero and 
one. Sin ce under EES the value of the exponent 7 only controls the time scale of creep iDanku and Kun 
(l2013al) the damage parameters were fixed to 7 = 1 and a = 1. 


When the load is put on the bundle, first some week fibers break immediately which may generate further 
breakings until a stable configuration is reached where all remaining intact fibers can sustain t he elevated 
load. The tim e evolution of the creeping system starts from this partially failed configuration [Kun et aT] 
(120071 [ 200 ^ . The present setup of the model impli es that the critical stress \ yhere immediate catastrophic 
failure occurs is equal to its static value ac = 0.25 iKun et alJ (120071120091120081) . In the presentation of 
the results the constant external load level will be characterized by the ratio as = (Tq/ ac which can take 
values in the range 0 < (Ts < 1 . 


2.2 BURSTS DRIVEN BY DAMAGE SEQUENCES 


The separation of time scales of the slow damage proce ss and of immediate bre a king leads to a highl y 
complex time evolution in agreement with experiments IKun et afl (I2009L120081) : iHalasz et alJ (120121) : 
Danku and KunI (l2013al]bh : starting from the initial configuration damaging fibers break slowly one-by- 
one gradually increasing the load on the remaining intact fibers. After a certain number of damage 
breakings the load increment becomes sufficient to induce the immediate breaking of a fiber which in turn 
triggers an entire burst of immediate breakings. As a consequence, the time evolution of creep rupture 
occurs as a series of bursts corresponding to the nucleation and propagation of cracks, separated by silent 
periods of slow damaging. The size of bursts A is defined as the number of fibers breaking in avalanches. 
The number of fibers breaking in a damage sequence and its duration determine the length of the damage 
sequence A^ and the physical waiting time T between consecutive events, respectively. 


Figure [T] shows a representative example of the time series of bursts as the system evolves towards 
failure. In order to have a clear view on the details of the sequence of events we intentionally used a 
relatively small system of A = 10® fibers subject to the load as = 0.001 which gave rise to 9793 
bursts. The size A of bursts is shown in the figure as a function of the discrete time, i.e. order number or 
natural time n = 1,2,... so that no information is presented about the physical waiting time T elapsed 
between the bursts. Strong fluctuations of the burst size A can be observed which is caused by the 
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Figure 1. Time series of bursts in a mean field (ELS) simulation of the FBM of creep. For demonstration purposes a relatively small system N = 10® is 
presented at the load as = 0.001 giving rise to 9793 avalanches. In order to have a clear view on the structure of the time series, sub-sequences starting at 
record breaking events are highlighted by different colors. The horizontal lines serve to demonstrate the increments between consecutive record bursts. 


quenched disorder of fiber strength. The time series has a non-stationary behavior which is indicated by the 
increasing average size of bursts. Of course, the physical waiting time T decreases between consecutive 
events, however, this info rmation is not visible in this representation . (The a nalysis of waiting times of 


the model can be found in iKun et alj (I2008L120091) : I Danku and KunI (I2013al) .) 


A record of the time series is a burst which has a size A^. larger than any previous events. Consecutive 
records are identified by their increasing rank /c as fc = 1, 2, 3,... which occurred as the n^th burst of the 
complete time series. The first burst n = 1 is by definition considered to be a record of rank k = 1 with 
ni = 1. In the example of Fig. [T] all together 18 record breaking events are identified (A: = 1,..., 18), 
which are highlighted by using different colors for the consecutive smaller events. It can be observed that 
RB events form a sub-sequence of bursts with monotonically increasing size, however, both the record 
size and the number of avalanches between two consecutive records exhibit strong fluctuations. In order 
to characterize these features we introduce the size increment 6k and the waiting time between two 
records with the definitions 


4 = A^+^ - A^, and mk = Uk+i - 


( 2 ) 


respectively. Note that the catastrophic burst which breaks all remaining intact fibers and destroys the 
bundle is not included in the time series so that the last burst in the bundle may not be an RB event. 


3 RESULTS 

To investigate the occurrence of record breaking events during the rupture process computer simulations 
were carried out for a system of size N = 10^ fibers averaging over 10^ realizations of the threshold 
disorders at each load value as- The external load as was varied over a broad range 0.001 < ag < 0.9, 
where the limits were set to have a sufficient number of bursts in the time series. We identify all record 
breaking avalanches A^ that occur up to macroscopic failure of the bundle and carry out a detailed analysis 
of their statistics. 
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Figure 2. (a) Fraction of fibers breaking in bursts {d^urst) due to damage {ddam) function of load Gs together with the number of bursts (A^a)- The 
quantity (dbuj-st) was calculated as the sum of the size of all bursts that occurred up to failure. The straight line has slope 0.87, furthermore, the vertical dashed 
line indicates the position of the characteristic load crj = 0.48. (b) Ratio of the average size of bursts (A) and of damage sequences (A^) as a function of the 
order number n of bursts for several load values, (c) Average size of the largest record breaking event and the average value of the largest waiting 

time between consecutive records as function of load. All quantities are normalized by the number of fibers N in the bundle. The upper and lower 

straight lines have slope 0.625 and 0.65, respectively. 


3.1 NUMBER OF RECORD BREAKING AVALANCHES 


An important feature of our system is that under low load values most of the fibers break due to slow 
damaging since the resulting load increments are too small to trigger extended bursts of immediate 
breakings. Figure[2](a) shows that the fraction of fibers breaking in avalanches (dburst) and due to damage 
{ddam) are monotonically increasing and decreasing functions of the external load ag, respectively. The 
two curves intersect each other at the characteristic load a* ~ 0.48, which coincides with the position 
of the maximum of the average number of avalanches {N^). The result demonstrates that for ag < 


cr: 


damage dominates the failure process, while the vicinity of catastrophic failure ag > a* is co ntrolled by 


the bu rsting activity with large burst sizes A and a decreasing number of bursts (see also iKun et al, 

dMl)). 

Figure [3](a) presents the average number of records (Nn) that occurred until n avalanches have been 
generated in the time series at several load values. Due to the breaking dynamics described above at low 
loads avalanches remain small typically comprising a few breaking fibers, and hence, new records mainly 
occur close to macroscopic failure. For instance, at loads ag < 0.01 the number of records has very 
low values (Nn) < 3 up to large n followed by a fast increase in the vicinity of the failure point. As 
ag increases the qualitative form of the (A„) curves remains the same, they just shift to higher record 
numbers due to the more intensive triggering of larger bursts at high loads. The most important feature 
of the record number {Nn) is that it has a logarithmic dependence on n over a broad range, which is in 
agreement with the analytic prediction of the logarithmic d epend ence of record numbers of IIDs on the 
event number n I Arnold et alJ (1201 Ih : [Shcherbakov et all (l2013h . The result implies that except for the 
close vicinity of macroscopic failure disorder dominates the process of creep rupture, and the occurrence 
of breaking avalanches can be well approximated as a stochastic process of IIDs. However, close to failure 
the increase of the stress on single fibers results in enhanced triggering which in turn gives rise to a sudden 
increase of records. The complete curves of (Nn) (n) can be characterized by the functional form 


(Nn) = A + B Inn + Cexp (n/D)^ 


(3) 


where the exponential term describes the rapid generation of high rank records close to catastrophic 
failure. It can be observed in Fig. [3][^a) that Eq. (|3]) provides an excellent fit of the numerical data over 
the complete load range considered. The additive parameter A has the same value A = 0.38 for all loads, 
while the multiplication factor B increases from 0.16 to 0.87 in the load range considered (see Fig. [3](6)). 
The shape of the {Nn) curves in the accelerating regime is controlled by the exponent ^ which is nearly 
constant ^ ^ 1.5 for loads below a* then it decreases to 0.33 as ag increases in Fig. [3]^5). The scale 
parameter D mainly sets the transition point between the two regimes dominated by disorder and by 
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Figure 3. (a) Average number of records (Nn) as a function of the total number of avalanches n for several load values. For low record numbers logarithmic 
dependence is evidenced, which is followed by a faster exponential increase. The continuous black lines represent fits with Eq. m- (6, c) Fitting parameters 
B, and D of Eq. as function of the external load as- (d) Probability distribution of the total number of records occurred up to failure for several 
load values. The distributions obtained at different loads are rescaled with the corresponding average and standard deviation. The continuous line represents 
the standard Gaussian Eq. 0. 


stress enhancement. It can be seen in Fig.[3lc) that D increases with the external load and has a maximum 
in the vicinity of a*, which is consistent with the behaviour of the total number of bursts (N^) in Fig. 
Efa) reflecting the overall dominance of damage and bursts in the failure process on the two sides of a*. 
To support the emergence of a characteristic event number D separating different regimes of bursting 
activity, Fig. ^b) presents the ratio of the average size of bursts (A) (n) and the average size of the 
damage sequence (A^) (n) that initiated the nth event. At the beginning of the process a large number 
of fibers must break randomly one-by-one in damage sequences until a small burst of size 1-2 fibers is 
initiated giving rise to (A) / (A^) < 1. Comparing to Fig.^a) the onset of the rapid exponential increase 
of record breakings corresponds to the event number from where extended bursts are triggered by shorter 
damage sequences due to the effect of stress enhancements. As the consequence of the increasing burst 
size and shortening damage sequences it becomes more and more probable that a record gets broken after 
a fewer bursts leading to the exponential increase of the record number. 

Figure [3l6) shows that the probability distribution ag) of the total number of records 

accumulated up to failure has a Gaussian form. The distributions ag) obtained at different loads 

were rescaled by the corresponding average and standard deviation cr^iot of the total number of 

records which results in a high quality collapse in the figure. The scaling function has a good agreement 
with the standard Gaussian 


p{x) 



exp [-x‘^/2), 


(4) 


which has been predicted for RB sequences of IIDs [Arnold et all (120 llh : [Shcherbakov et alJ (120131) . The 
Gaussian functional form prevails in spite of the complex effect of the acceleration of the occurrence of 


Frontiers in Interdisciplinary Physics 


6 




























Danku et al. 


Record breaking in creep rupture 




Figure 4. Size distribution of record breaking bursts (a) and their increments (c) for several load values, (b) and (d) demonstrate that rescaling the 
distributions with the external load high quality collapse of the curves can be achieved. The continuous lines in (6) and (d) represent fits with Eq. 0. 


records towards failure. Deviations can be observed at the tails of the distribution which may indicate a 
slight right-handed asymmetry. 


3.2 STATISTICS OF RECORD SIZES AND WAITING TIMES 


Recently, it has been shown that in our model the size distribution of avalanches v(A ) accumulating all 
crack ling events up to failure during the creep process follows a power law distribution iKun et alJ (12008 . 


200 ^ 


p(A) ~ A 


(5) 


where the value of the exponent r coincides with the usual mean field exponent of FBMs r = 5/2 
Kloster et alJ (Il997l) . As t he external load a pproaches the critical value cts —)■ 1 a crossover is obtained to 


a lower exponent r = 3 /2 IKun et alJ (l2008h in agreement with the mean field prediction ofsimple FBMs 


subject to a quasi-statically increasing load lPradhan et ^ (1201 OL l2005h : iRaischel et alj (120061) . 


Figure Sla) presents the size distribution p(Aj., as) of record breaking bursts accumulating all records 
up to failure at different load values. The distributions can be well described by a power law followed by 
a finite size cutoff of exponential form. Since at higher loads larger avalanches are triggered the cutoff 
of the distributions shifts to higher values but the functional form remains the same for all loads. Figure 
|4f6) illustrates that rescaling the two axis in Fig. Ufa) with appropriate powers of the external load a high 
quality data collapse can be achieved. The result implies that record size distributions obtained at different 
loads obey the scaling structure 


p{Ar,as) = A^^’’(/(Ar/cT"), 


where the exponent a controls the load dependence of the cutoff record 


(Amax) 


a 


a 

s • 


( 6 ) 

(7) 


Frontiers in Interdisciplinary Physics 


7 
















































Danku et al. 


Record breaking in creep rupture 



Figure 5. (a) Probability distribution of waiting times p(m, as) between consecutive records for several values of the external load as- In [b) and (c) the 
same distributions are presented after rescaling the two axis with powers of the external load as separately above and below a*, respectively. The good quality 
data collapse demonstrates the validity of the scaling structure Eq. m- 


The power law dependenee of the eutoff reeord size (AJ?®®) is also eonfirmed by Fig. Efc), where the 
largest reeords of single simulations were direetly averaged. Best fit is obtained with the exponent 
a = 0.65 whieh provides also the best eollapse in Fig. Note that due to the normalization of the 
distributions the sealing exponent used for the transformation along the vertical axis has to be equal 
to the product of Tr and a which is explicitly indicated in Fig. ^h) . The size distribution exponent of 
record bursts = 1.33 ± 0.0 3 proved to be significantly s malle r than the usual niean field bu rst size 
exponents of FBMs iKun et al] (I2008L l2009t) : iKloster et alJ (Il997h : Pradhan et alj (I2005L l2010ll . which 
is the consequence of the relatively high frequency of large burst sizes in the RB sequence with respect 
to the complete time series. The short flat region for the smallest bursts occurs due to down sampling of 
small sized events in the RB sequence. To fit the scaling function in Fig.|4](6) the cutoff function 0(a;) was 
assumed to have a stretched exponential form 4>{x) ~ exp(—( x/xq)^) with (3 = 0.65. It is important to 
emphasize that the distribution p{Ar, as) has a homogeneous evolution with increasing external load, i.e. 
the exponent Tr remains constant when the critical load is approached ds —)■ 1. 


For the advancement of the RB sequence of bursts during the evolution of the rupture process the size 
increments 6k = — A^ between consecutive records carry also interesting information. Figure IH^c) 

presents that the distribution p(<5,., as) of increments 5r has a qualitatively similar behavior to record sizes 
p{Ar, as), i.e. power laws are obtained followed by an exponential cutoff as described by Eq. ®. Careful 
scaling analysis in Fig. Wid) shows that both the value of the exponent Tr of the power law regime and 
the scaling exponent a of as have the same values for the two distributions, the only difference is that the 
cutoffs of increments are smaller than the one of record sizes. 


After a record occurred as the n^th avalanche of the time series it gets broken after a certain number 
of events by the njt+ith avalanche. The waiting time between RB avalanches defined by Eq. Q is 
an important characteristic quantity of the RB sequence of bursts. Eigure|5l^a) presents the accumulated 
statistics of waiting times m considering all re cords k in the RB sequence for several load values. Based 
on the statistics of extremes Galambo^ (119781), for IIDs a power law behavior is expected p{rn) ~ m~^ 
with the exponent zjib = 1 lArnold et alJ (l201 Ih . It can be observed in the figure that our results are 
generally consistent with the IID prediction, however, distributions p(m, as) obtained below and above 
form two groups of different power law exponents. Eigures[5|&) and (c) demonstrate that rescaling m 


a 


and the distributions p{m, as) with appropriate powers of the external load as, good quality data collapse 
can be achieved in both load regimes. The scaling structure is similar to Eq. 


p{m) = m ^ijj{m/ag), (8) 

however, both exponents 2 ; and a proved to have different values below and above the characteristic load 
a*. In Eigures [5](5) and (c) best collapse is obtained with the exponents z = 0.72, a = —1.45 (a* < as) 


Frontiers in Interdisciplinary Physics 


8 




































Danku et al. 


Record breaking in creep rupture 


and z = 1.15, a = 0.625 (a* > ds). The value of z refleets an interesting aspeet of the dynamies of the 
ereep proeess: the low value of z < zjjd at high loads shows that long waiting times more frequently 
oeeur than for IID. The reason is that due to the large bursts triggered under high external loads it takes 
longer for the system to break a reeord. At low loads the waiting time exponent 2 ; is slightly larger than 
the IID predietion implying an elevated frequeney of short waiting times with respeet to IIDs. Note that 
a has different signs in the two regimes eorresponding to the inereasing and deereasing behavior of the 
eutoff of the distribution with inereasiim load. An independent test of the load dependenee of the eutoff 
waiting time is shown in Fig.|2|c), where was direetly averaged over the simulations. A 

maximum is obtained at a* in agreement with the sealing behavior in Figs. [5]^ 6) and (c). The flattening of 
the distributions p(m, ag) for low m values and the small bump elose to the eutoff observed for low loads 
in Fig.[5]^c) are eaused by the finite system size and by the distinet distribution of the time needed to break 
the first reeord, respeetively. 


3.3 EVOLUTION OF THE SEQUENCE OF RECORDS 


The ereeping system approaehes maeroseopie failure through an aeeelerating sequenee of bursts of 
inereasing size. In order to understand how reeord breaking events oeeur during this evolution we 
evaluated average values of the eharaeteristie quantities of single reeords as a funetion of their rank k. 
In Fig.[^a) the average reeord size (A^) has the same generie form for all loads, i.e. a nearly exponential 
inerease is obtained with a slight minimum of the derivative of the eurves for intermediate ranks. Note that 
for a given value of /c a reeord ean have a higher value at low loads than at higher ones, e.g. at fc = 15 the 
eorresponding reeord burst is signifieantly larger at (jg = 0.001 than at cxs = 0.1 in spite of the two orders 
of magnitude higher external load in the seeond ease. The reason is that for avalanehe triggering the load 
on single fibers is the most relevant quantity, henee, at the generation of the fcth reeord the ereeping system 
ean be eloser to eatastrophie failure at lower loads than at higher ones, whieh implies larger burst sizes. In 
order to gain information about the rate of inerease of reeords, we determined the relative inerements of 
eonseeutive RB events defined as the ratio 5fc/A^. Figure shows that at the beginning of the rupture 
proeess the average (^^/A^) starts from a high value simply beeause the first reeord has size Aj = 1 
and it gets typieally broken by a burst of size A^ = 2 or A^ = 3. Then the relative inerement rapidly 
deereases to the vieinity of 0.25. The remarkable result is that for reeords of the highest rank (5^/ 
tends to 1, whieh implies that as the system approaehes maeroseopie failure reeord breakings oeeur by 
nearly doubling the size of the previous reeord. Henee, not only the reeord sizes but also the inerements 
form a monotonieally inereasing sequenee when approaehing failure. 


As the exter nal load inereases larger a valanehes are triggered in the system and they oeeur with an 
inereasing rate iPanku and Kunl (12013al) . Based on this general tendeney it eould be expeeted that the 
average time (rifc) where the k\h reeord appears is a deereasing funetion of the external load ag. Figure 
[^c) shows that for low load values just the opposite happens: reeords of a given rank k oeeur later after a 
larger number of avalanehes. The reason is that as the load inereases early reeords get larger, and henee, 
they are more diffleult to overeome. The situation ehanges at the eharaeteristie load a* so that in the range 
a > a* the average reeord time (rifc) deereases with inereasing load. 


The time series presented in Fig.[T]does not eontain information about the physieal time elapsed between 
events. However, the aeeeleration of the system ean still be inferred from the natural time n, beeause 
the number of bursts between RB events deereases when approaehing eatastrophie failure. Due to this 
aeeeleration the average waiting time (mjt) is expeeted to deerease between eonseeutive reeords. It is 
interesting to note that for higher loads the deereasing braneh of (mfc) is preeeded by a rapidly inereasing 
regime indieating the slow down of reeord breaking before aeeeleration sets on. The maximum of {m^ 
already develops for the lowest loads with the position k* = 2 and it gradually shifts to k* = 13 for the 
highest ones. Comparing Figs.|^6) and (d) it ean be seen that the reeord rank k* of the maximum of (rrik) 
falls elose to the position of the minimum of the relative inerements (^6k/^r)- The result demonstrates the 
emergenee of a eharaeteristie time seale rik* of the system whieh separates slow down and aeeeleration of 
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Figure 6. Characteristic quantities of single records: the average value of the size (Aj:) (a), and relative increments of record breaking events (b). 

In (b) the horizontal dashed line highlights the limit value 1. Time of occurrence (rife) (c) and waiting time (mfe) between consecutive records (d). 


the dynamics of record breaking. This is also in agreement with the behavior of record numbers in Fig. 
H a) so that n^* approximately coincides with the onset time of the exponential increase of {Nn). 


4 DISCUSSION 

We investigated the statistics of records in a sequence of crackling avalanches which occur during the 
creep rupture of heterogeneous materials. Synthetic sequences of bursts were generated by computer 
simulations of a realistic fiber bundle model where slowly developing damage triggers avalanches of 
immediate breaking of fibers. The bundle is subject to a constant external load below the fracture strength 
of the system. We analyzed the mean field limit of the model where all fibers keep the same load so that 
no spatial correlation develops between local failure events. 

Record events are defined solely based on the burst size, i.e. a record is a bursts whose size is larger 
than that of all previous bursts. This way a monotonically increasing sub-sequence of crackling events 
is identified. Computer simulations revealed that during the evolution of the rupture process the average 
number of records increases logarithmically with the number of avalanches except for the close vicinity 
of macroscopic failure where an exponential form is evidenced. Additionally, the total number of records 
obtained up to failure at different load values has a Gaussian distribution. These findings are in agreement 
with the robust analytic predictions on the RB statistics of sequences of independent identically distributed 
random variables which shows that the beginning of the creep process is mainly controlled by the 
quenched disorder of the system. The enhanced triggering close to failure due to the rapidly increasing 
load on single fibers is responsible for the exponential acceleration of record numbers. 

The size of records proved to have a power law distribution with an exponent 1.33 significantly lower 
than the usual mean field exponents 5/2 and 3/2 of the burst sizes of FBMs. The size increments between 
consecutive records are found to have the same scaling structure and the value of exponents as the record 
size. To prove the independence of the exponents on the external load a careful data collapse analysis 
was performed. The probability distribution of waiting times between consecutive events has also a power 
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law functional form, however, with different exponents below and above the eharaeteristie load a*. In 
order to eharaeterize the evolution of the sequenee of reeords we studied the average value of the relative 
inerement and of the waiting time between eonseeutive events as a funetion of the reeord rank. Both 
quantities show the emergenee of a load dependent eharaeteristie seale in the system: at the beginning a 
slow down of reeord breaking oeeurs with an inereasing waiting time and deereasing relative inerement. 
Beyond a eharaeteristie reeord rank k* the approaeh to failure results in an aeeeleration of reeord breaking 
with deereasing waiting times and inereasing relative inerements. In our study the disorder of materials 
was represented by uniformly distributed failure thresholds. We repeated the eomplete RB analysis for 
Weibull distributed disorder varying the Weibull exponent. It has to be emphasized that the qualitative 
behavior of the results and the value of the exponents of the reeord size, inerement and waiting time 
distributions all proved to be universal, only the eharaeteristie load a*, furthermore, the sealing and eutoff 
exponents a and [3 depend on the disorder. 


Reeently, the reeord breaking statis ties of driven threshold models of eomplex systems has been 


analyzed in IShcherbakov et alj (120131) . These are eellular automata models of self organized eritieality 
where the slow external driving leads to the emergenee of a steady state eharaeterized by intermittent 
avalanehes of relaxation events. Comparing the statisties of reeords of the models to the eorresponding 
results of IIDs the authors eould point out eorrelations in the eomplex spatio-temporal evolution of 
avalanehes. The mo st prominent deviation from IIDs was found for the Olami-Feder-Christensen model 


Olami et alJ (119921) where the number of reeords proved to inerease as a power of the logarithm of the 


avalanehe number. The main differenee of the dynamies of our FBM and the above models is that during 
ereep both the size and the rate of avalanehes inerease so that no steady state arises. Henee, the eomparison 
to IIDs in our ease helped to determine regimes eontrolled either by disorder or by the inereasing stress 
level of intaet fibers. Reeord b reaking statis t ies of inter-event times in aftershoek sequenees of earthquakes 
has reeently been studied in [Yoder et^D (120101) . Based on the non-homogeneous Poissonian proeess 
deseribing the rate of events, a power law behavior of the reeord number was obtained following the 
logarithmie inerease. Our study demonstrates that the investigation of the reeord breaking statisties of the 
time series of eraekling events reveals also interesting novel aspeets of rupture phenomena. 
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